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ABSTRACT 

We present a preliminary analysis of the small-scale structure found in new 70-520 yum continuum maps of the Rosette molecular cloud (RMC), 
obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure 
identification algorithm, with sizes up to ~ 1 .0 pc in diameter. A comparison with recent Spitzer maps reveals that 37 1 clumps are "starless" (without 
an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have 
masses (Mc) over the range -0.75 < log (McfM^) < 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes 
that are consistent with those found for high-mass clumps identified in CO emission by other groups. 
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1. Introduction 

Stars form within molecular clouds, after some fraction of cloud 
material is first condensed into smaller-scale structures. How this 
process unfolds is not well understood, however. For example, 
what conditions within clouds drive the formation of small-scale 
structures that in turn produce stars of any given mass, e.g., high- 
mass stars. Detailed observations of the small-scale structure 
within molecular clouds should provide valuable insight into 
how stars of various masses form out of dense material. 

The need for sensitive probes of small-scale structure in 
molecular clouds is the impetus behind the "Herschel OB Young 
Stellar objects" (HOBYS) Key Project (see Motte et al. 2010). 
The HOBYS team is currently using the ESA Herschel Space 
Observatory (Pilbratt et al. 2010) to obtain wide-field maps of 
~15 molecular clouds which are forming high-mass stars within 
3 kpc of the Sun. These clouds are being observed with both 
the Herschel SPIRE (Griffin et al. 2010) and PACS (Poglitsch et 
al. 2010) instruments, ultimately to obtain wide maps of these 
clouds at 70-520 /im at diffraction-limited resolutions of ~18" x 
(A/250 /im). Such data sample thermal emission from cold dust 
grains mixed with molecular gas within the cloud. Since contin- 
uum emission at submillimetre wavelengths is typically optically 
thin, detections of continuum emission from dust at high resolu- 
tion can trace the organization of mass in clouds on small scales. 
By sensitively sampling emission across the peaks of the spectral 



Send offprint requests to: J. Di Francesco 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
Correspondence to: james.difrancesco@nrc-cnrc.gc.ca 



energy distributions (SEDs) of this dust, SPIRE and PACS data 
can constrain both the column density and temperature of the 
dust and provide unparalleled censuses of the small-scale struc- 
ture in molecular clouds. 

In this paper, we describe the small-scale structure detected 
with SPIRE and PACS in the Rosette molecular cloud (RMC), 
the first of the HOBYS target list that was observed as part of 
the early Herschel "science demonstration phase" (SDP) cam- 
paign. The RMC is a rich location of star formation at a dis- 
tance of 1.6 kpc (Johnson 1962; Perez et al. 1987; Park & Sung 
2002). The cloud is southeast of NGC 2244, the Rosette Nebula, 
and the expanding HII region associated with NGC 2244 has 
begun to interact with it. The structure of material within the 
RMC has been previously examined by Williams, Blitz & Stark 
(1995), Schneider et al. (1998) and Dent et al. (2009) using var- 
ious CO observations. Alternatively, the populations of young 
stellar objects (YSOs) within the RMC have been studied by 
Phelps & Lada (1998), Li & Smith (2005), Poulton et al. (2008), 
and Roman-Zuiiiga et al. (2008) using infrared observations. 
Notably, these studies identified 10 embedded clusters within the 
RMC, i.e., PLOl-07 and REFL08-10. A recent review of obser- 
vations of the RMC has been made by Roman-Zuiiiga & Lada 
(2008; RZL). 

The SPIRE and PACS observations described in this paper 
provide a preliminary look into how Herschel data trace small- 
scale structures in the RMC, and the relationships between these 
structures and YSOs. Analyses of other aspects within the SDP 
data of the RMC can be found in the accompanying papers by 
Motte et al. 2010, Henneman et al. 2010, and Schneider et al. 
2010.) 
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2. Observations 

A complete description of the observations described here can 
be found in the accompanying paper by Motte et al. (2010). 
In short, an extinction map of the RMC was first constructed 
for planning purposes using 2MASS data and the AvMAP algo- 
rithm. AvMAP derives extinction maps with typical resolutions 
of ~2' using colour excesses in a manner akin to the NICER al- 
gorithm of Lombardi & Alves (2001). AvMAP, however, also 
uses a stellar model to predict foreground source densities; see 
Schneider et al. (2010a). Using the AvMAP data as a guide, an 
area of the RMC ~ 1 deg^ in extent encompassing areas south- 
east of NGC 2244 with Ay > 5 was mapped in common with 
SPIRE and PACS in parallel scan mode on 2009 October 20. 
Maps were obtained simultaneously at 70 /im and 170 /vm with 
PACS and at 250 yum, 350 yum, and 520 fj.m with SPIRE. 

The SPIRE and PACS data were reduced separately us- 
ing common routines of the Herschel Interactive Processing 
Environmen|J(HIPE; see Ott 2010). To facilitate the identifica- 
tion and analysis of the cold, small-scale structure in the RMC, 
i.e., dense clumps, the PACS maps were Gaussian smoothed to 
the 18" FWHM resolution of the 250 ^iva SPIRE map. A full 
analysis of the unsmoothed PACS maps will be provided in the 
near future by Sadavoy et al. (2010b, in preparation). 

3. Results 

Figure 1 (available electronically) shows the 350 yum map of the 
RMC, revealing a rich amount of structure in continuum emis- 
sion. (See Motte et al. for maps at all five wavelengths.) Several 
bright knots are seen embedded in large filaments, and each can 
be associated with known infrared clusters (see Figure 1, also 
RZL). Ubiquitous, faint diffuse emission is also seen, and in 
many cases this appears filamentary as well. Of special note in 
Figure 1 is the filamentary arc seen in the northwestern corner, 
indicating the "ridge" of the RMC {+ "Core D"), i.e., where the 
adjacent NGC 2244 HII region (the Rosette Nebula) is interact- 
ing with the molecular cloud. 

From the (smoothed) PACS and SPIRE maps, 473 contin- 
uum objects were identified using the getsources algorithm; 
an accompanying paper in this volume by Men'shchikov et al. 
(2010) describes getsources in more detail. In brief, getsources 
first decomposes maps into images over a range of scales and 
objects are identified as discrete emission regions in each im- 
age. Objects are then defined using information on many scales 
and at many wavelengths. For example, positions and extents 
(i.e., "footprints") are determined from the first and second 
moments of intensities in images where all scales are recom- 
bined. Integrated flux densities are calculated by summing the 
background-subtracted emission under the footprint of an object. 
Sizes are well-approximated by 2D Gaussian fits. 

For simplicity, we choose to identify all objects found by 
getsources as "clumps'nFigure 1 shows the locations of all 473 
clumps identified in the RMC overlaid onto the 350 //m map. 
Most are found within the brighter filaments of the RMC, with 
neighbouring clumps tracing the linear structures of underlying 
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^ The term "clump" has been used previously to describe sites of 
grouped or clustered star formation, as opposed to sites of single- (or 
multiple-) star formation, i.e., "cores"; see Williams, Blitz & McKee 
(2000). We use "clump" here to describe all discrete condensations of 
material of size < 1 pc within the larger cloud. 



filaments. Several "isolated" clumps (i.e., those without close 
neighbours), however, are also seen associated with fainter emis- 
sion. Relatively few clumps are seen inside the northwestern 
"ridge" and coincident with the NGC 2244 HII region. 

Figure 1 also shows which clumps in the RMC contain YSOs 
("protostellar clumps") and which do not ("starless clumps"). 
These determinations were made using Spitzer IRAC and MIPS 
data of the region (see Poulton et al. 2008; Balog et al. 2010, in 
preparation). The criteria used to determine which clumps con- 
tain YSOs are based on those developed recently by Sadavoy 
et al. (2010a) for small-scale structure detected at 850 yum in 
the less-distant Orion, Taurus, Perseus and Ophiuchus molecular 
clouds. For this comparison, potential YSOs had to be detected 
in all four IRAC bands or at 24 jum. Their colours and brightness 
were checked against colour criteria to distinguish them as ei- 
ther YSOs or other objects (e.g., background galaxies or stars). 
In addition, their positions relative to clumps were checked for 
coincidence within areas defined by the 2D Gaussian fits to the 
clump extents. We do not distinguish here between protostellar 
Classes of YSOs in clumps; see Henneman et al. for further dis- 
cussion. 
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Fig. 2. Distributions of 350 /vm fluxes (A) and geometric mean 
sizes (B) of RMC clumps. Dotted black, solid red, and solid 
blue lines indicate all clumps, starless clumps, and protostel- 
lar clumps, respectively. In panel A, the vertical dashed line in- 
dicates roughly a 3 o" rms limit for detected clumps. In panel 
B, the vertical dashed line indicates the resolution limit of the 
smoothed 250 yum data. 



In total, we find 371 starless clumps and 102 protostellar 
clumps in the RMC. Most protostellar clumps contain one YSO, 
but ~35% contain more than one. Figure 1 suggests that brighter 
clumps tend to be protostellar, while fainter clumps tend to be 
starless. (Exceptions to both cases are seen, however.) As ex- 
pected, the protostellar clumps are associated with the embed- 
ded clusters known in the cloud, i.e., PLOI-07 and REFL08-10. 
Many starless clumps appear to be associated with the RMC 
ridge and the area between it and the RMC "central core" (i.e., 
the concentration of bright protostellar clumps associated with 
the embedded clusters PL04-06 and REFL08 in the centre of 
Figure 1). Relatively few (~15) clumps are seen behind the ridge 
itself, and they all appear to be starless. Excluding areas inside 
and outside the ridge, the starless and protostellar clumps ap- 
pear well-mixed, suggesting no "age gradient" across the RMC. 
(Schneider et al. 2010b suggest, however, that the most massive 
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cores in the RMC decrease in age with increased distance from 
NGC 2244.) 

Figure 2 shows distributions of the flux and size of clumps 
in the RMC, as a whole and when separated into starless or pro- 
tostellar clumps. Figure 2a shows histograms of 350 /im inte- 
grated fluxes, 5 350, of the clumps, demonstrating how the bright- 
est are indeed protostellar (see Figure 1). Starless and protostel- 
lar clumps with S 350 ~ 5 Jy are similar in number, but clumps 
fainter than 5 Jy are predominantly starless. For reference (see 
§4 below), a clump of 5 Jy flux at 350 //m in the RMC has a 
mass of ~7 Mq, assuming Tdusi = 20 K. Figure 2b shows size 
histograms of each population where size is the geometric mean 
of the observed major and minor axes of the 2D Gaussian fit to 
each clump extent. (These sizes are not deconvolved from the 
finite beam.) Clumps of each type have sizes that extend from 
the resolution limit up to ~0.65 pc. The few larger clumps (i.e., 
0.7-1.0 pc) are all protostellar, however. 
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Fig. 3. Distributions of local visual extinctions (A) and tempera- 
tures (B) for the RMC clumps. Lines defined as for Figure 2. 

Figure 3 shows distributions of the local extinction and tem- 
peratures of clumps within the cloud, again as a whole and when 
separated into starless or protostellar clumps. Figure 3a shows 
the Ay values associated with the clumps, as determined from 
the low-resolution (i.e., ~2') extinction maps. For a given clump, 
its associated low-resolution Ay can be used as a proxy for its 
depth within its parent cloud (see Johnstone, Di Francesco & 
Kirk 2004). Clumps of each type are found across the same range 
of extinction {Ay ~ 3-20). No clumps of either type are found at 
Ay < 2, however. Starless clumps are found in greater numbers 
at lower extinctions (Ay < 12), but both types are found in simi- 
lar numbers at higher extinctions (Ay > 12). The number of star- 
less clumps peaks at Ay ^ 6, while that of protostellar clumps 
peaks at Ay x 13-14. Though starless and protostellar clumps 
are found over the same range of Ay, the difference in peak Ay 
may be due in part to protostellar clumps having higher masses 
on average than starless clumps (see §4 below). Assuming sim- 
ilar sizes, the resulting relatively higher column densities could 
mean relatively higher local Ay values, though the diff'erence in 
the AvMAP and Herschel resolutions is large. 

Figure 3b shows the distributions of the clump temperatures 
(Tdust) found across the RMC using all five Herschel wave- 
lengths and a method described by Schneider et al. (2010b) 
Temperatures could be found for only 402 of the 473 identified 
clumps (i.e., 323 starless and 79 protostellar clumps), and we 



limit further discussion only to those clumps. Dust temperatures 
for these clumps range between 10 K and 30 K. Interestingly, 
starless clumps appear to be generally warmer than protostellar 
clumps, a surprise as the latter could be expected to be warmer 
given internal heating by YSOs. The difference here, however, 
stems from the large population of starless clumps located to the 
northwest of the RMC (see Figure 1). At this location, the dust 
is quite warm owing to its proximity to the O stars in NGC 2244. 



4. Clump mass spectra 

The mass for each clump, Mc, was determined using its ob- 
served 350 /urn flux and the relation Mc = S yd" I KyByiJ du^i) 
= 0.073(5 v/Jy)(^/kpc)2(^,/0.07 cm2g-i)-i[exp(41.1/rrf„,,)-l] 
Mq. To account for variations in temperature across the RMC, 
we assume for each clump its respective Tdust value when deter- 
mining its mass, rather than averages of the distributions shown 
in Figure 3b. We also use kt,s() - 0.07 cm^ g^'for a dust opacity, 
consistent with the dust prescription of Ossenkopf & Henning 
(1994), and a distance of 1.6 kpc. 

Figure 4 shows the distributions of mass for the total, star- 
less, and protostellar clump populations, i.e., their clump mass 
spectra (CMS). Masses for all clumps range across ~3.5 orders 
of magnitude, i.e., -0.75 < log (Mc/Mq) < 2.50. Evidence for dif- 
ferences in mass between the populations is seen, however, with 
starless clumps having on average lower masses than protostellar 
clumps. Moreover, the starless population peaks at log (Mc/Mq) 

- 0.25, while the protostellar population peaks at log (Mc/Mq) 

- 1.0. This diff'erence is likely because the starless clumps in the 
RMC are both warmer and have lower submillimetre fluxes on 
average than the protostellar clumps (see Figures 2 and 3). Note 
also that there are fewer higher mass starless clumps than proto- 
stellar clumps, e.g., over the high-mass range of log (Mc/Mq) = 
1 .75-2.50, only two starless clumps are found but 1 1 protostellar 
clumps. 
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Fig. 4. Clump mass spectra. Lines defined as for Figure 2. 
For comparison, the solid black line indicates the slope of the 
steeper, high-mass tail of the Salpeter IMF. 

Table 1 lists the linear slopes (a) fit to the high-mass tails 
of the respective CMS, using the least-squares method and as- 
suming Poissonian errors in each mass bin. Table 1 also lists the 
mass ranges used to compute the slopes. In general, slopes were 
fit between the bin containing the peak of the respective CMS 
and the bin of the highest mass in the sample. The a - -0.65 ± 
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0.01 for all high-mass clumps is slightly more than the value of 
-0.8 found by Dent et al. (2009) from CO 3-2 observations over a 
5 deg^ map of the entire Rosette complex. The value found here, 
however, is quite similar to the -0.6 slope found by Schneider 
et al. (1998) for higher clump masses using CO (2-1) maps of 
the Rosette. Both these groups used different clump identifica- 
tion algorithms than used here, i.e., clumpfind by Dent et al. and 
gaussclumps by Schneider et al. (1998). Moreover, we have de- 
termined masses from clumps we identified here, not over the 
extents of the clumps found by either group using CO. 



of 10-30 K, with starless clumps surprisingly warmer on average 
than the protostellar clumps, owing to the proximity of many to 
NGC 2244. Masses were determined for each clump, revealing 
a range of -0.75 < log (Mc/M©) < 2.50. Linear least-squares fits 
to the high-mass tails of the CMS reveal slopes that are consis- 
tent with slopes fit to mass spectra of clumps identified through 
CO in the RMC. Such slopes are shallower than those seen in 
the high-mass tails of the mass spectra of small-scale continuum 
structures ("cores") in closer clouds, possibly from crowding of 
such objects on scales smaller than probed here. 



Table 1. Slopes of high-mass tails of RMC clump mass spectra 



Clump Type 



Mass Range 



All -0.65 ± 0.01 0.25 < log (Mc/Mq) < 2.50 

Starless -0.82 ± 0.02 0.25 < log (Mc/Mq) < 2.25 

Protostellar -0.80 ±0.02 1 .00 < log (Mc/Mq) > 2.50 



After dividing the clumps by type, our starless and protostel- 
lar clumps have similar slopes to their high-mass distributions, 
i.e., a - -0.8, but over different mass ranges. Indeed, these slopes 
are more similar to that found by Dent et al. The high-mass tails 
of the all-clump CMS has a shallower slope than do just starless 
or protostellar clumps. This difference stems from the difference 
in the mass ranges used to measure slopes. Also, the numbers 
of protostellar clumps peak at a higher mass than do those of 
starless clumps, and they dominate the all-clump CMS at even 
higher masses. 

With single fits, all three CMS are shallower at higher masses 
than seen in the high-mass tails of mass spectra of small-scale 
structures within closer clouds. Indeed, the slopes of such tails 
are seen to be consistent with the -1.35 slope of the Salpeter IMF 
(e.g., in Aquila; see Andre et al. 2010). For comparison with the 
RMC, a Salpeter slope is also shown in Figure 4. This differ- 
ence may come simply from resolution; namely, in closer clouds, 
sites of individual star formation, ("cores") can be resolved. For 
example, typical core separations in Ophiuchus at 125 pc were 
found to be -0.03 pc by Motte, Andre & Neri (1998). Though 
our minimum linear resolution (at 250 jum) is ~0. 12 pc, or about 
the expected size of isolated cores, we may not be able to resolve 
individual cores in the RMC owing to crowding on scales < 0.12 
pc. Indeed, such crowding may cause the resulting mass spectra 
to be artificially shallow. Higher resolution data will be critical 
for testing this possibility. 

Without appropriate spectral line data, we cannot estimate 
the states of gravitaitonal stability for the clumps in the RMC, 
and determine which may be bound or which may be transient 
structures. For example, we cannot presently guess which star- 
less clumps will likely engage in future star formation. It is clear 
from Figure 1, however, that the Herschel data have revealed that 
star formation is ongoing in the RMC and that we are seeing only 
a snapshot of a very dynamical situation. 

5. Summary 

We have conducted a preliminary analysis of small-scale emis- 
sion within the Rosette molecular cloud, as observed by the 
Herschel Space Observatory. From SPIRE and smoothed PACS 
images, we identified 473 clumps using getsources, a new al- 
gorithm that identifies structure over multiple wavelengths and 
scales. Using Spitzer data, we classified these into 371 starless 
and 102 protostellar clumps. The clumps have dust temperatures 
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Fig. 1. Herschel map of the Rosette molecular cloud at 350 /im. The positions of 473 clumps identified by getsources are shown, 
where red crosses denote starless clumps and green boxes denote protostellar clumps. The cyan dashed line indicates the approximate 
edge of the Rosette Nebula and the cyan cross indicates the central position of NGC 2244 from Table 3 of RZL. The nine embedded 
clusters identified by Phelps & Lada (1998; PL) and Roman-Ziiniga et al. (2008; REFL) in the Herschel field are also labelled. 



